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Results from ultra-relativistic Au+Au collisions at RHIC are reviewed. Emphasis is put upon (a) measured 
properties of a Au+Au collision, which might be used as input to cosmic air shower Monte-Carlo event generators, 
(b) production of anti- matter, and (c) forward physics. 



1. RHIC 

The Relativistic Heavy Ion Collider (RHIC) at 
Brookhavcn National Laboratory, New York, 
USA, has a circumference of 3833 m and uses 
1740 superconducting magnets pQ. Fig. ^ 
shows the RHIC accelerator complex. The 
main objective of RHIC is the investigation 
of novel QCD phenomena at high density and 
high temperature in ultra-relativistic Au+Au 
collisions. The highest center-of-mass energy 2 
^=200 GeV is a factor ~10 higher than past 
fixed target nucleus- nucleus collision experiments 
(e.g. CERN SPS E beam =W0 GeV, equivalent 
to ^^17.3 GeV). With a cross section of 
^Au+Au—7-04: barn (see Tab. 2), the design lu- 
minosity of £=2-10 26 cm _2 s _1 means an interac- 
tion rate of ~1.4 kHz. Four experiments are lo- 
cated at RHIC interaction points, i.e. the STAR, 
PHENIX, PHOBOS and BRAHMS experiments. 
The STAR experiment consists of 

• a large scale midrapidity (|?y|<1.6) Time 
Projection Chamber (TPC, R=2 m, L=4 m) 
with ~48,000,000 pixels (momentum resolution 
Apt/pt=3% at Pt=1 GeV/c in a solenodial mag- 
netic field of B=0.5 T) [2] 

• 2400 lead-scintillator electromagnetic calorime- 
ter modules (|?7|<1.0), and 

• a 3-layer silicon drift detector with ~13,000,000 
pixels (M<1.0) 
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Figure 1. The RHIC accelerator complex. 

The PHENIX experiment 3 consists of 

• a midrapidity (|ij|=±0.35) electron and hadron 
detection system, using a ring imaging Cerenkov 
detector, a time-of-flight system, drift chambers 
(mass resolution Am/ra=0.4% at m=l GeV) and 
24,768 lead-scintillator electromagnetic calorime- 
ter modules, and 

• forward/backward (1.1<|?/|<2.2) muon cham- 
bers. 

The PHOBOS experiment g] uses 135,168 sili- 
con strip and silicon pad readout channels in a 
coverage of 1 77 j ==1=5-4, A<^=360°, for unidentified 
charged particles, and \r)\=±2.0, A^=180°±20° 
for identified charged particles. 
The BRAHMS experiment is a spectrometer 
with a small opening angle of AO=7.3 msr, how- 
ever at an adjustable rapidity position, thus pro- 
viding particle identification from |?y|=0 to 1 77 1 =4 . 
As an example for all experiments, Tab. 1 shows 
the number of recorded events at the STAR ex- 



1 



2 



periment since start of RHIC operation. 



Table 1 

Recorded STAR data set. 







(J of Events 


Year 


Au+Au 


130.0 GeV 


0.7 Mill. 


2000 


Au+Au 


200.0 GeV 


3.2 Mill. 


2001 


Au+Au 


19.6 GeV 


-20.000 


2001 


d+Au 


200.0 GeV 


35 Mill. 


2003 


Au+Au 


200.0 GeV 


77.9 Mill. 


2004 


Au+Au 


62.4 GeV 


13.3 Mill. 


2004 


p+p 


200 GeV 


19.8 Mill. 


2004 



Au+Au collisions create a system of hot nuclear 
matter, d+Au collisions a system of cold nuclear 
matter. The p+p data set serves as reference for 
Au+Au, but at the same time is also used for the 
RHIC spin program, as protons in RHIC can be 
polarized vertical or longitudinal. 
As an illustration of the amount of data, the raw 
data set collected at STAR experiment in 2004 
compares to ~60 times the combined raw data set 
of the B meson factories BELLE and BABAR. 
STAR also utilizes a level-3 trigger system [B], 
performing full event reconstruction of a cen- 
tral Au+Au collision within i<100 ms. Fig. 
shows a central Au+Au collision at -^=200 GeV, 
recorded by the level-3 trigger system at STAR, 
and consisting of ~6,500 charged particle tracks 
and -130,000 TPC clusters. The level-3 trigger 
was used for realtime rare probe detection, such 
as events with anti-Helium (described in detail in 
Ch.EJ 

2. The RHIC Au+Au Collision 

There are four basic stages of an Au+Au collision, 
which are also schematically depicted in Fig. 3. 

• t=0 fm/c. The nuclei are Lorentz contracted 
in the laboratory frame ( "pancake shape" ) . Hard 
parton scattering occurs (qq, qg, gg scatterings). 

• i~l fm/c. A hot cylinder is formed. The tem- 
perature in the cylinder exceeds 10 12 K. 

• i~4 fm/c. Soft partons, thermally produced, 
emerge from the center of the collision. Beam 
remnant particles are peaked in the forward and 
backward directions. 

• £~10 fm/c. Hadronic freeze-out occurs. All 



partons confine into particles. 




Figure 2. RHIC Au+Au collision at 

Vs=200 GeV, recorded at the STAR exper- 
iment. 

Table 2 

Measured parameters of a RHIC Au+Au collision. 
For details see Ch. [21 









Rcf. 


GAu+Au 


(7.05±0.05) barn 


130 




N part 


352±3 


200 


8 


Nbin 


1051±72 


200 


8 


Rout 


(5.39±0.18±0.28) fin 


130 


9 


-Rlong 


(5.99±0.19±0.36) fm 


130 


9 




(5.48±0.15±0.30) fm 


130 


9 




(620±33) GeV 


200 


A 


Tpein 

LL/rp 


(216±14) GeV 


200 


m 


P 


0.56c 


130 


ma 


T^chera 


(181±8) MeV 


130 


na 




(89±10) MeV 


200 


EH 


Mb/T 


0.18±0.03 


130 




p/p ratio 


0.71±0.002±0.05 


130 


d 


£ Bjorken 


(4.2±0.3) GeV/fm 3 


130 


na 



Tab. 2 lists the measured parameters of a RHIC 
collision. pau+Au denotes the total cross sec- 
tion. N part denotes the number of participant 
nuclcons. Nti n denotes the number of binary 
nucleoli- nuclcon collisions. The size of the sys- 
tem can be concluded from the radii R ou t, Riong 
and Rside, which are measured by quantum- 
mechanical Hanbury-Brown-Twiss (HBT) inter- 
ferometry of 7r ± 7r ± pairs, close in phase space. 
Rout denotes the radius perpendicular to the 
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beam axis, Ri ong parallel to the beam axis. R S ide 
denotes a radius perpendicular to R ou t, but not 
necessarily measured from the origin (0,0,0). R out 
can be interpreted as the size of the particle wave 
front, while freezing out. In particular, the ratio 
Rout/Rside— 1 indicates that RHIC freeze-out has 
a surprisingly short duration, an evidence which 
is supported by the measured high expansion ve- 
locity, given as the fraction of the speed of light 
,3=0.56. 

Et denotes the total transverse energy, defined as 
the sum of the electromagnetic transverse energy 
Erf 71 (i.e. leptons and photons) and the hadronic 
transverse energy E T ad . The ratio E T m /Et was 
determined to be 0.348±0.019 0. So far, there 
is no indication for Centauro type events with an 
anomalous ratio E^ n /E T ad ~l/Q fTS] . 
The p/p ratio visualizes how many baryons from 
the beam projectiles remain in the collision zone. 
A complete baryon-free zone would lead to a 
p/p ratio ~1, i.e. all baryons and anti-baryons 
are created by pair production. The measured 
p/p ratio of 0.71 indicates that a significant ex- 
cess of baryons over anti-baryons is still present 
at RHIC. As a comparison, the measured ra- 
tio is much higher than 0.00025±10% at the 
AGS (v/s-4.9 GeV) and 0.07±10% at the SPS 
(•y/i~17.3 GeV), i.e. RHIC collisions represent for 
the first time an opportunity to advance into the 
almost baryon free regime. 
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Figure 3. The four basic stages of an ultra- 
relativistic Au+Au collision. For details see 
Ch.H 

The baryo-chemical potential /ib denotes the 



energy which is necessary to add one nucleon 
to the system. Recent lattice QCD calcula- 
tions indicate the position of the tri-critical point 
of strongly interacting matter at a temperature 
T=(160±3.5) MeV and a baryo-chemical poten- 
tial of /i b =(725±35) MeV [TB]. Thus, the primor- 
dial conditions in a RHIC collision are close to the 
critical temperature, but the baryo-chemical po- 
tential has a difference of a few hundred MeV to 
the tri-critical point 3 . As one consequence, large 
scale event-by-event fluctuations in the number of 
charged hadrons (as one might expect close to the 
tri-critical point) are not expected at RHIC. How- 
ever, it should be noted, that in the 2-dimensional 
QCD phase diagram (T vs. /i&) the primordial 
conditions of a RHIC collision are very close to 
those of the primordial universe. 



2.1. The Temperature 

One distinguishes between the chemical freeze- 
out temperature T c h em and the kinetic freeze-out 
temperature Tki n , which correspond to different 
times in the system evolution. t(T c h em ) denotes 
the end of the inelastic collisions inside the par- 
ticle wave front, i.e. relative particle abundances 
become fixed. t(Tj~i n ) denotes the end of the elas- 
tic collisions and is ~100 MeV cooler than T c h em - 
The temperature T c fc, em =181 MeV corresponds to 
2.1-10 12 K, which can be compared to the temper- 
ature in other systems (Tab. 3). The temperature 
resembles one of the highest temperatures in the 
universe. 

Also an interesting fact from the compari- 
son of the temperatures at y / s=130 GeV and 
•^5=200 GeV shall be noted: although the col- 
lision energy increases by AE—70 GeV, the tem- 
perature only increases by AT— 7 MeV, which 
corresponds to a small fraction of <10~ 6 A- AE . 



3 Note that the new GSI accelerator facility FAIR (see 
Ch.|HJ will test a (/ij,,T) region in the vicinity of the pre- 
dicted tri-critical point, where event-by-event fluctuations 
might occur. 
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Table 3 

Comparison of temperatures in various systems. 

1.4-10 34 K Planck temperature 

2.1-10 12 K Au+Au collision 

~10 n K mini black hole (R=l fm, M=10 12 kg) 

~10 9 K supernova 

15- 10 6 K sun (core) 

55- 10 6 K plasma fusion 

4-10 6 K laser fusion 

3T0 4 K thunderstorm flash 

2.2. The Matter Density 

The spatial matter density, i.e. the number of par- 
tons per volume p=N part ons/V can be estimated 
by a Bjorken ansatz |17j . If assuming that ini- 
tial and final entropy are equal, the number of 
partons at t<l fm/c is equal to the measured 
number of final state hadrons Nhadron- The vo- 
lume can be calculated by inserting the fireball 
radius {R ou t and Ri ong in Tab. 2) into a cylin- 
drical volume (due to Lorentz boost in the beam 
direction), i.e. p~dN part on/ dy-1 / (nR 2 t) using the 
rapidity y. For a time early in the system evolu- 
tion i=0.2 fm/c the matter density is j o~20/fm 3 , 
which corresponds to ~15xp D , p denoting the 
density of cold gold nuclei. At this high density, 
hadrons are definitely non-existent. In Tab. 4, 
the density is compared to other systems. 

Table 4 

Comparison of the density in various systems. 



30-10 17 kg/cm 3 
2-10 17 kg/cm 3 
~20,000 kg/cm 3 
-1000 kg/cm 3 

l.l-lO" 26 kg/cm 3 



Au+Au collision 
Au nuclear density 
Au atomic density (solid) 
metallic hydrogen 
(core of planet Jupiter) 
universe critical density 



2.3. The Energy Density 

About 90% of all emerging particle are 7r mesons, 
thus one may assume (as an educated guess) the 
7r meson mass in the particle energy E 2 =m 2 +p 2 . 
By (a) counting the number of charged parti- 



4 The rapidity y of a particle is defined as 
y=l/2-\n((E+p z )/(E+p z )) with the total particle 
energy E and the momentum component in beam 
direction p z . The yield dN/dy is Lorentz invariant. 



cles, (b) measuring the average particle momen- 
tum p and (c) using the volume according to 
Ch. 12.21 one can estimate the total energy den- 
sity as e~5 GeV/fm 3 , corresponding to a factor 
~30 higher energy density then in cold Au nu- 
clei. For comparison, the estimate for the AGS 
is e=1.2 GeV/fm 3 , for the SPS e=2.4 GeV/fm 3 
|18j . A better estimate based upon perturbative 
QCD can be found elsewhere [7], leading to an 
even higher energy density of e~18 GeV/fm 3 . 

3. From p+p to Au+Au 

As RHIC experiments measure both Au+Au 
and p+p collisions in the same detector, cor- 
responding data of both collision types can be 
directly compared with identical systematic er- 
rors. Fig. 0] shows transverse momentum distri- 
butions for 7r~, K~ and p for ^=200 GeV p+p 
(top) and Au+Au (bottom) as a compilation of 
all four RHIC experiments. One can immedi- 
ately see that a Au+Au collision is not a sim- 
ple superposition of many p+p collisions. In fact, 
the average mean transverse momentum <pt> is 
~390 MeV/c for p+p and ~508 MeV/c for cen- 
tral Au+Au. The increase by ~30% indicates the 
presence of collective effects. For further details 
see Ql] [H]. 

4. Qualitative comparison to the primor- 
dial universe 

A RHIC Au+Au collision is sometimes being ref- 
ered to as the "little bang", in reference to the 
"big bang" of the primordial universe. However, 
a few important differences shall be noted: 

• The primordial universe underwent a matter- 
dominated and a radiation-dominated period. At 
RHIC, matter dominates exclusively. 

• At the time when the universe reached RHIC 
temperatures of T~100 MeV, it already had 
developed a macroscopical horizon distance of 
L~10 km, and thus had a size which was roughly 
larger than the size of a RHIC collision by a factor 
10 18 . 

• The expansion velocity of w=0.56c at RHIC in- 
dicates an almost explosive character. In case 
of the primordial universe (after inflation) it was 
much slower by a factor 1 / Mpi anc k~10 19 . 
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"beam" energies) range from 10 1 -10 6 MeV per 
nucleon, comparable to RHIC beam energies at 
10 5 MeV per nucleon. A particularly suited probe 
for anti-matter search is anti-He, as the signature 
of charge Z=—2 is clean 5 . A Z——1 signature of 
p is contaminated by K~ and 7r~ mesons. The 
STAR level-3 trigger system was used to identify 
Z=—2 candidate events in realtime. Primordial 
Helium is not rare, but in fact represents ~24% of 
all primordially created matter (created at time 
scales W s and at temperatures T~0.1 MeV). 
Under the assumption of an a priori symmetric 
primordial production of He and anti-iJe, one 
may conclude that anti-ffe might be detectable, 
too. 

In a RHIC Au+Au collision, there is trivially 
no anti-matter in the initial state, as nuclei are 
completely made of matter. However, even anti- 
nuclei are created in a detectable quantity. The 
dominant anti-nuclei production mechanism is a 
two-step process, namely (1) pair production of 
pp, nn, followed by (2) coalescence, i.e. the anti- 
nucleon wave functions overlap inside a homogen- 
ity volume. For the systematic study of the A 
dependence of anti-nuclei production yields, it is 
useful to define an invariant yield 



Figure 4. Transverse momentum distribution for 
7r~, K~ and p for • v /s=200 GeV p+p (top) and 
Au+Au (bottom). For details see Ch.|3 

• In the system evolution of the universe, RHIC 
matter density of p~15p Q occured only for a 
very short period during inflation (assuming in- 
flation from 10 +118 p o at £~1CT 36 s to KT 19 ^ 
at i^l0~ 34 s). Additional details are given else- 
where |21| . 

5. Anti-Helium 

Currently one of the most interesting questions 
of cosmology is: where is the primordial anti- 
matter ? Several experiments (e.g. the bal- 
loon bound BESS experiment [22] at a height 
of H=27 km and the space-shuttle bound AMS 
experiment at an orbit of H— 320-390 km) 
search for direct evidence of anti-matter in cos- 
mic rays. It should be noted that incident parti- 
cle energies at the AMS orbit location (i.e. quasi- 
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with anti- nuclei yield N Al the anti- nucleon yield 
Nn, the anti- nucleon momentum p, and p=P / A. 
The coalescence coefficient Ba denotes the prob- 
ability that A anti-nucleons form a bound state. 
It might be regarded as a "penalty factor" for the 
step from n to n + 1 anti-nucleons in the system. 
Typical orders of magnitude are i?2~10~ 3 and 
-B3^10 -6 . For high y/s and a large system size, 
the coalescence coefficient is related to the inverse 
of the effective volume containing the anti-nuclei 
(i.e. the fireball volume) by Ba—1/V^'J 1 . 
So far, STAR has been able to collect a raw yield 
of 193 anti- 3 i/e in a total of 4-10 6 Au+Au events 
at -^5=200 GeV. In the same data set, 6416 anti- 
deuterons were also identified and used for a de- 



5 Except anti-He and hypothetical pentaquark states there 
are no other stable Z=—2 particles predicted in nature. 
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termination of B 2 . Further details are given else- 
where [2U- Fig. 5 shows the measured ionization 
dE/dx vs. momentum for the anti--ffe candidates 
in the STAR Time Projection Chamber. 
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Figure 5. Ionization dE/dx vs. momentum p for 
anti-ffe in the STAR Time Projection Chamber 
(blue : level-3 trigger, red : offline analysis). 

Fig. 6 shows the extracted coalescence coefficient 
£?3 in comparison to other measurements. RHIC 
represents the highest energy ever at which anti- 
He was created in a accelerator experiment. 
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Figure 6. Coalescence coefficient B3 for anti- 3 £fe 
and anti-t production vs. y/s for different experi- 
ments. 

As an important result, the coalescence (and thus 
the yield of produced anti-nuclei per one collision) 
in Au+Au (normalized to the number of binary 
collisions) at RHIC is approximately a factor of 
40 smaller then in p+p or in 7+p 25 . 



The coalescence coefficient is also an impor- 
tant input to calculations of anti-TJe yield from 
nucleus-nucleus collisions in the stellar medium 
|26| . which is the dominant background for the 
primordial anti-ffe search by AMS. 
The STAR 2004 data set is currently being ana- 
lyzed. The data are expected to contain 2.7±0.3 
anti- 4 7Je, which would be a first experimental ob- 
servation. 

6. Forward Physics 

For high energy p+A collisions in the earth's up- 
per atmosphere, a proper knowledge of the parti- 
cle shower in the forward rapidity region is very 
important, as the forward particles are the visible 
center of the showers as seen on the detector sys- 
tems at sea level. Ideally p+Au collisions would 
be preferable, but due to technical accelerator is- 
sues 6 d+Au collisions were chosen at RHIC. As 
shown in Tab. 1, several million d+Au collisions 
were recorded at RHIC in 2003. Fig. 7 shows the 
rapidity distribution of particles from a d+Au col- 
lision at ^5=200 GeV with data from BRAHMS 
|27j . in comparison to a perturbative QCD cal- 
culation ■ The 77 asymmetry ratio of forward 
(Au side) to backward (d side) is a characteristic 
of the collision system. STAR measured the ra- 
tio for these collisions as a function of pr (Fig. 8) 
|29j. As a surprising result, for 0< 1 77 1 <0 . 5 the ra- 
tio is flat, but for 0.5<|?7|<1.0 the ration shows 
a rising behaviour up to pt— 2-3 GeV/c, followed 
by a descending slope for higher pf. Fig. 8 also 
shows the comparison to model calculations, in- 
cluding (a) shadowing 7 and (b) saturation at high 
gluon densities. The data seem to support a satu- 
ration ansatz, which, in a geometrical picture, as- 
sumes that the probability that two gluons collide 
is one. p+A collisions will also be an important 
part of the future GSI accelerator facility FAIR 
at Darmstadt, Germany (Fig. 9), which will pro- 
vide p beams for fixed target experiments up to 



6 One the one hand, the d beam has a charge/mass ratio 
close to the Au beam. On the other hand the beam-beam 
interaction, which shortens the lifetime of a RHIC store, 
is in the case of p+Au a factor of ~2 higher than in the 
case of d+Au. 

7 Shadowing is the depletion of parton distributions at 
small x inside a nucleus. 
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Figure 7. Pseudorapidity distribution for charged 
particles from d+Au collisions at v / s=200 GeV 
in BRAHMS ;27 for two different centralities, 
corresponding to two different impact parameter 
ranges. 
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Figure 8. Forward/backward pseudorapidity ra- 
tio vs. pt d+Au collision at -^5=200 GeV in 
STAR For details see Ch.® 

Another important class of forward reactions are 
ultra-peripheral Au+Au collisions, in which the 
Au nuclei interact by a very strong electromag- 
netic field (10 20 V/cm at the surface of the nu- 
clei). Photons produced in these strong fields are 
peaked forward, and may create vector mesons 
by scattering off the Au nuclei. Fig. 10 shows 
the differential cross section d<r( r yAu^p°Au)/dt 
vs. the momentum transfer t. The very small 
t values indicate the forward scattering. From 
an exponential fit (shown as line in Fig. 10) 
a forward cross section was determined to be 




Figure 9. 
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Figure 10. Forward differential cross section 
da('yAu->p Au)/dt vs. t [23]. 

dcr/ffi| t=0 =965±140±230 mb/GeV 2 . Details are 
given elsewhere |30| . 

7. Charm cross sections 

For high energy cosmic nuclear collisions in the 
galactic medium or the earth's atmosphere, the 
knowledge of heavy quark production cross sec- 
tions is very important. Both STAR jJH] and 
PHENIX J>2] have measured the charm cross sec- 
tion. Fig. 11 shows the total cc cross section per 
nucleon-nucleon collision vs. y/s. There is a pre- 
liminary indication that the measured cross sec- 
tion is higher than the expected cross section ex- 
trapolation obtained from (a) PYTHIA Monte- 
Carlo simulation (based upon Lund fragmenta- 
tion) and (b)& NLO QCD calculation 
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Figure 11. Total cc cross section per nucleon- 
nucleon collision vs. y/s. F° r details see |31| . 

8. Summary 

In 4 years of operation, RHIC has collected a 
data set of millions of p+p, d+Au and Au+Au 
collisions. From the analysis we may conclude 
that a RHIC Au+Au collision is a very hot and 
very dense hadronic system, which expands ex- 
plosively with v=0.56c. During the short freeze- 
out the system undergoes a graceful QCD self- 
organisation, which could be successfully de- 
scribed by models based upon perturbative QCD. 
One of the most basic but also most important 
observations is that a Au+Au collisions is not 
a simple superposition of p+p collisions, but in 
fact the mean transverse momentum increases by 
~30%, indicating collective effects. 
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